Mapping interacting onto non-interacting quantum Hall systems 
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We construct an explicit duality between the interacting quantum Hall system in the lowest 
Landau level and a non-interacting Landau problem. This is done by absorbing the interaction into 
the gauge field in the form of an effective magnetic vector potential. The result is analogous to, 
and illuminates the microscopic origin of, the well-known composite fermion model, but has several 
advantageous properties. Using this duality we derive, for an arbitrary short-range interaction, 
analytic expressions for the ground state energy and the excitation gap as functions of the filling 
fraction. We find good agreement with existing results. 
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Our theoretical understanding of the fractional quan- 
tum hall effect is based on a number of simple, but very 
successful qualitative descriptions. The most prominent 
of these are the Laughlin trial wave-functions [l| and 
the composite fermion picture of Jain 0] ■ Despite these 
successes there is, apart from the work of Murthy and 
Shankar j3] , no true microscopic derivation of these de- 
scriptions. Here we attempt to construct such a dual 
description explicitly. 

We consider a system of N spinless electrons moving in 
a two dimensional region with surface area A. The parti- 
cles experience an uniform magnetic field perpendicular 
to the plane and interact via a repulsive two-body poten- 
tial V{r). The Hamiltonian reads 



H 



-y 



(1) 



where e = — |e| is the charge and m the mass of an elec- 
tron. Interparticle distances are denoted by = jr^ | = 
|ri — rj|. In the symmetric gauge the vector potential 
for an uniform magnetic field Bz is — {B/2)z x r^. 
It is convenient to treat V {rij ) as dimensionless and 
let Vo have the dimension of energy. The kinetic en- 
ergy term in _ff is a regular Landau problem which can 
be solved exactly Q and produces the well-known pic- 
ture of degenerate Landau levels separated by a gap of 
tujjc where ojc — {\e\B) / {mc). Let 0o = he/ \e\ and 
a = \J (he) I (|e|_B) denote the flux quantum and mag- 
netic length respectively. The degeneracy of each Landau 
level equals AB and the filling fraction v = N(j>o/AB 
is the, typically fractional, number of levels filled by the 
electrons. In what follows will always denote the filling 
fraction of the original Hamiltonian in Although it 
does not appear explicitly in H, we assume that a neu- 
tralizing background charge is present in the system. 

Following the literature we restrict ourselves to the 
lowest Landau level (LLL). When H is projected onto 
this subspace the kinetic energy is a trivial constant, leav- 
ing only the projection of the interaction: 



(2) 



Later we will encounter Landau problems where singular 
flux attachments [H| have been made to the particles. 
We use Hq {B + J |) to denote such a Hamiltonian with 
an uniform magnetic field Bz and flux attachments that 
carry 7 flux quanta each. 

Now we introduce a new Hamiltonian H, which is a 
modified Landau problem with an additional term added 
to the vector potential. The form of H is 



H 



(3) 



where Ai = z x ^E*"* (r;) and 



E*°*(r,) = ^E(r,,)=^£;(r,,)f.,. (4) 

Here 7 = {'^(j>o) / {27Ta) and E (r) is a scalar function. 
Note that both E*°* (r) and the parameter 7 are dimen- 
sionless. Let us consider the magnetic field corresponding 
to Ai ■ The magnetic field at due to z x 7E (r^ ) is 

B (r„-) = z • [V X (z X 7E (ry))] = 7V • E (r„-) , (5) 



while the total field at generated by Ai is 



5*°* (r,) = 7V • E* 



(6) 



We interpret this by saying that Ai attaches a magnetic 
field to each particle. The radial profile of this field is 
determined by E{r), e.g. when E{r) = 1/r these are the 
familiar singular flux attachments [5|. 

We will show that VHV is closely related to the pro- 
jected interaction in ([2]) when Ai is defined according to 



E{r)r^ -aVV (r) 



(7) 



For short-range interactions E{r) will decay quickly at 
large r. This implies that the total flux associated with 
Ai is zero, since by Stokes' theorem J B • da ~ E{r)r and 
E{r)r — ^ as r — > cxD. When the particles are well sep- 
arated Ai will therefore have no effect. This reflects the 
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short-range nature of the interactions which Ai was con- 
structed to mimic. 

In order to perform some necessary simphfications wc 
will use a mean-field approximation based on an uniform 
particle density j5 = N/A = v/{2'Ka?'). Due to the inter- 
actions and fermion statistics the particles are strongly 
correlated, and we incorporate this in our approximation 
using the following simple picture. Imagine that around 
each particle there is a region of radius in which the 
probability of finding another particle is very small, taken 
in our case to be zero. The notion that each particle oc- 
cupies a certain minimum area is supported by the fact 
that within the LLL the projected spacial coordinates of 
a particle no longer commute This is implemented 
in the mean-field approximation by taking the particle 
density, as seen by a specific particle at r^, to be 



|r - TjI < re 
p otherwise 



(8) 



i.e. a uniform density beginning at a distance re from r^. 
The excluded region around is denoted by A\. To fix 
re we argue that at low energy the system would tend 
to maximise the interparticle separation. This effectively 
assigns an area of A/N — I-ko^ jv to each particle. Tak- 
ing this region to be a disk of radius we conclude that 



a 



(9) 



We remark that the precise value of re is not crucial, and 
does not affect the duality's qualitative results. 

One application of this scheme is to approximate the 
attached fields by a single uniform field. For physical 
interactions V(r) should decay smoothly at r > re. We 
expect B(r) to share this property, which lends itself to 
a mean-field approximation: 



B'°\v,) = Y^B(r,,) ^ -pi drB(|r,-r|) 



(10) 



where A = reE{re)/a is an interaction dependent param- 
eter. The integral over A was dropped since B{r) carries 
no net flux. The flnal expression follows from applying 
Stokes' theorem to the integral over A*. 



Now we turn to the main task of projecting H onto 
the LLL. Let V denote the projection operator onto the 
subspace spanned by Slater-determinants formed using 
the LLL states of Hq{B). Take VHV and multiply out 
the Ai term to obtain 



4i„- + 4f 



(11a) 

r (lib) 



+72^PE*°*(rO •E*°*(r,)7' (11c) 

i 

where Eq = hujc/2 and = {2h)^^z ■ [r^ x pij] is 
the z component of the dimensionless relative angular 
momentum operator. Since Hq(B) assumes a constant 
value of NEq within its LLL the first term (|llap is 
simply N. The remaining two terms are treated sepa- 
rately, starting with (jllb|l . Here we use the fact that 
any operator containing only relative coordinates is com- 
pletely characterized in the LLL by its matrix elements 
between states of the same relative angular momentum. 
These are the Haldane pseudopotentials f7] and the rele- 
vant states are \n) — {zi — zj)"' exp [— rj^/4a^ — rj^/4a^] 
with zj. — Xk — iyig. and Lij \n) = —n\n). Note that 
fermion statistics require n to be an odd positive in- 
teger. It is straightforward to check that (|llb[) has 
the same pseudopotentials as, and is therefore equiva- 
lent to, (27ra2/(/)o) J2i VB*°^ [vi) V. We wiU approximate 
Qtot ^^^^ j-|y mean- field value derived in (jlOp . This 
reduces (|llbp to a simple scalar: —v^XN. 

Next consider line (jllcp and rewrite the sum as 



j,k L i 



r,: - r,- 



E(r, -rfe) 



(12) 

We will use the mean-field approximation to replace the 
sum over i by an integral. Following ([8]) the regions 
around and are excluded from the domain. The 
expression above then becomes 



A-Ai-Ai; 



drE(r-rj)-E(r-rfe). 



(13) 
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Next we rewrite the dot product using the identity E (r — rj) • E (r — r^) 



-V • [aV (r - rj)E(r - r*.,) 



aV (r — Tj) V • E (r — Tfe) and then apply Stokes' theorem to the integrals of total divergences. This leads to 



y [E*°* (r,)] » -ap i da ■ (r) E*°* (r)l -I- — / drl/*°* (r) (r) 

, JdA 190 J A 



E 



apf da- [F(r -rj)E(r -rfe)] +ap i da • (r - r^) E (r - r^)] 

JdAi JdAi 



(14) 
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Here V*°* (r) and (r) are the total potential and 
magnetic field at r due to particles which are at least 
a distance of from r. We will consider each term in 
Pl)) separately. The integral over OA may be neglected 
based on the following electrostatic analogy. Equations 
([7]) and ([5|) suggest that we may consider E*°* (r) as the 
electric field due to a charge density distribution propor- 
tional to (r). Since 5*°* (r) carries no net flux the 
total charge in this analogy is zero. The integral over 
dA therefore involves the electric field at the edge of a 
rotationally invariant system with zero net charge. By 
Gauss' law E*"* (r) can then be taken as zero. Correc- 
tions that arise due to boundary effects will not scale 
extensively and is neglible in the thermodynamic limit. 
The integral over dA^ is proportional to V{re), which 
we can choose as zero through a trivial shift in V(r). 
This is also why no j = fc terms appear in the dou- 
ble sum in p^ . For the integral over dA^ we approx- 
imate V {r — Ti) by expanding it around r = as 
V{r-Vj) = V (r,fe) + E (r,,.) ■ (r - r^) /a. The sec- 
ond term of the expansion does not contribute to the 
integral, leaving lyXV {rjf^ ) . Finally consider the integral 
over A. Since i3'°' (r) only contains contributions from 
particles lying at least Ve from r it has the same mean- 
field value as B*°^ (r^) in (fTUI) . even though r is not a 
particle coordinate. The mean-field value of y*°*(r) is 
Vbg (r) — 2Ve where T4 — t^P Jq° drrV{r) is a constant 
and Vbg (r) — pj^ dr' V (r' — r) has the form of a po- 
tential due to a neutralizing background charge. We now 
expand the argument of the integral over A to linear order 
in the fluctuations around these mean- field values. The 
resulting expression can be treated using approximations 
very similar to those already encountered. Combining 
the results of this discussion leads to 



(15) 



where Hi„t contains the particle-particle, particle- 
background and background-background interactions: 

H^nt = ^ ^ inj) - J2 ^"9 (r.) + f / Vtg (r) . (16) 

Returning to (jlip and replacing each line by its simplified 
form gives 



V 



H 

Eq 



a'N 



r = 2pXj^P [H,nt + NVe] V (17) 



where a' — 1 — v^X. This is the desired relationship 
between VHV and the projected interaction. To exploit 
this result we need to solve VHV, which we do using a 
mean-field approach. 

First we turn H into a free particle Landau problem by 
replacing the attached fields by an effective uniform field 
with strength —vjEX, as shown in (fTO|) . This renormal- 
izes the existing field Bz by a factor of a' < 1 . The flux in 



the excluded region Al around ri due to the field attached 
to particle i cannot be probed directly by other parti- 
cles. It does however affect the Aharonov-Bohm phases of 
these particles. It follows that this flux acts as an effective 
flux attachment carrying A7 flux quanta. The mean-field 
approximation of H now becomes Hq {a' B + 7' ']) where 
a' — 1 — v^X and 7' — 7A. The relation a' = 1 — v^' en- 
sures that the total flux present in H and Hq {a'B + 7' "f ) 
are the equal. Upon replacing H by Hq {a'B + 7' t) 
equation ()17p becomes strongly reminiscent of the com- 
posite fermion picture 0, @| in that it relates a Landau 
problem with a weakened magnetic field plus flux attach- 
ments to an interacting quantum Hall system. For even 
integer values of 7' the Hamiltonian Hq {a'B + 7' |) is 
exactly solvable, although its projection is generally not. 
However, if the low energy states of Hq {a'B + 7' t) lie 
predominantly within the LLL of Hq {B), we may drop 
the projection operators on the left of P?|) . We argue 
that this is indeed the case. This can be understood in- 
tuitively by observing that if the average flux density in 
Hq (B) and Hq {a'B + 7' t) are equal, states which are 
localised on a radius r will enclose, on average, the same 
amount of flux and therefore have the same Arahanov- 
Bohm phase. This will maximise the overlap between 
the two sets of states. Our construction already sat- 
isfies the requirements of this argument. First, it fol- 
lows from a' — 1 — v^' that the total flux in Hq {B) and 
Hq {a'B + 7' t) are equal. Secondly, the uniform particle 
density ensures that the average density of singular flux 
attachments are also uniform, implying that the systems 
have an equal flux density. It follows that for low energy 
states with a uniform density, the projection does not 
dramatically affect the spectrum. The gap in particular 
should be preserved. We drop the projection operators 
on the left of p?|) and obtain the final form of the duality: 



Hq {a'B + y T) 

Eq 



-a'N = 



X 



V[H,nt+NVe]V. (18) 



The free parameter 7 was originally introduced in |3|. 
Its renormalized value 7' — A7 will be chosen as in the 
composite fermion literature, where the main considera- 
tion is the removal of ground-state degeneracies fsj. For 
the famous series of filling fractions i' — {2k + 1/p)"^ 
we set 7' — 2k, in which case the filling fraction of 
Ho{a'B + j' I) is lyQ = 1/p. The flux attachments do 
not alter the statistics of the particles, and the ground 
state corresponds to filling the lowest p Landau levels. 
The finite gap at this filling reveals the incompressibility 
of the quantum Hall ground state. Expressions for the 
ground state energy Eg and excitation gap A of VHintV 
can easily be read off from (fT8|) as 



E. 



a 
N 



-Ve + 



(19) 
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1/3 2/5 3/7 4/9 5/11 

Ai/^ ([21} 0.10206 0.0559 0.0385 0.0295 0.0238 
MDD [9] 0.1012 0.05 0.035 0.027 - 
CFT [10| 0.1005 0.0549 0.0371 0.0276 0.0219 



5/11 



6/13 



7/15 



Ai/^ ((22| 0.0238 0.02 0.0173 
SLJ fll] 0.0219(30) 0.0225(41) 0.018(11) 



TABLE I: The gaps for the Coulomb interaction from ([22] 
and BEm. 



the extrapolation method used in PMJ; an issue which 
was raised in MDD. Not shown in the tables are results 
for the Coulomb interaction at = 1/5 which appear in 
11 [S [13 and range from 0.023 to 0.025. Equation 
gives 0.0198; a difference of about 20%. 



V 


1/3 


2/5 


3/7 


1/5 


2/9 


Eg IN JH]) 


-0.4082 


-0.4193 


-0.4243 


-0.3162 


-0.3281 


JK [15] 


-0.4098 


-0.4328 


-0.4422 


-0.3275 


-0.3427 



TABLE IL The ground state energy 
Coulomb interaction from (|2ip and [IE 



per particle for the 



and 



A 



4Pp " (y-vf 



(20) 



We now identify the origin of the FQHE's robustness with 
respect to changes in the form of Y{r): By absorbing the 
interactions into the gauge potential in the form of at- 
tached fields a mean-field treatment is possible in which 
all information about the precise form of Vij) is lost. 
The interaction dependency only enters through A which 
rescales the spectrum but cannot affect the existence of 
the gap. To test the accuracy of and ([^ . we con- 
sider some specific examples. The ground state energy 
per particle for the Coulomb interaction (a/r) is given 
by 



N 



- 1 



(21) 



while the excitation gaps for the a^'-fr^ and 
aexp [r/(aK)] /r (Yukawa) interactions are 



A 



l/r" 



nv p f v 



id 



i/2 



Ay = 



a ^ 1 



(22) 
(23) 



Tables |T] to IIIII compare our analytic results to those ob- 
tained through various other schemes. These include 
numerical diagonalization as in MDD and the com- 
posite fermion theory (CFT) developed in [l^, [13, [I3|- 
PMJ Q, JK and SLJ [uj made use of Monte Carlo 
methods. The numerical values quoted for PMJ are esti- 
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mates based on the graphs in 

For the Coulomb interaction we found good agreement 
with the results of CFT, MDD and JK with discrepancies 
of less then 10% for the gaps and 5% for the ground state 
energy. Our results for the gap at 5/11, 6/13 and 7/15 
agree particularly well with those of SLJ. At these filling 
fractions up to seven composite fermion Landau levels 
are filled. At some filling fractions the results of PMJ 
differ significantly from ours for the 1/r^ and Yukawa 
interactions. This may be due to finite size effects or 



1/3 



2/5 



3/7 



4/9 



/r2 



22]) 0.0833 



0.05 



0.0357 



0.0278 



PMJ [14] 0.0842(11) 0.0609(27) 0.0424(45) 0.0257(77) 



1/3 



2/5 



3/7 



4/9 



Ay 1231 0.0304 0.01933 0.01406 0.011 
PMJ [14] 0.0271(4) 0.0215(9) 0.0143(15) 0.0102(24) 

TABLE III: The gaps for the 1/r^ (top) and Yukawa (bottom) 
interactions from ((^ . ([^ and 14]. We set k = 1. 
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